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^ ■ Abstract 

D 

Q^l A measurement is presented of elastic Deeply Virtual Compton Scattering e+ + p — ;• 

^ ' e"*" + 7 + p at HERA using data taken with the HI detector. The cross section is measured 

as a function of the photon virtuality, Q^, and the invariant mass, W, of the jp system, in 
the kinematic range 2 < Q^ < 20GeV^ 30 < 1^ < 120 GeV and \t\ < 1 GeV^ where 
/\ ' t is the squared momentum transfer to the proton. The measurement is compared to QCD 

j^ ■ based calculations. 
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1 Introduction 



This paper presents the first measurement of the elastic cross section for Deeply Virtual 
Compton Scattering (DVCS) (Fig. [I|a) i.e. the diffractive scattering of a virtual photon off a 
proton [jl]-^ by studying the reaction: 



P 



e^ + 7 + p. 



(1) 



The interest of this process, because of its apparent simplicity, resides in the particular insight 
it gives for the applicability of perturbative Quantum Chromo Dynamics (QCD) in the field of 
diffractive interactions. The wide kinematic range in the photon virtuality, Q^, accessible at 
HERA, provides a powerful probe for the interplay between perturbative and non-perturbative 
regimes in QCD. Furthermore the DVCS process gives access to a new class of parton distri- 
bution functions, the skewed parton distributions (SPD) [[7]-^ which are generalisations of the 
familiar parton distributions and include parton momentum correlations. 

The reaction studied receives contributions from both the DVCS process whose origin lies 
in the strong interaction, and the purely electromagnetic Bethe-Heitler (BH) process (Fig. [T]b 
and c) where the photon is emitted from the positron line. The BH process is precisely known as 
it depends only on QED calculations and proton elastic form factors. The DVCS cross section 
is extracted by subtracting the BH contribution from the total cross section, which is possible 
since the interference contribution vanishes when averaged over the full azimuthal angle of the 
final state particles. A recent measurement of the single spin asymmetry in a longitudinally 
polarised electron beam [ ]T0[ ] complements this measurement. 
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Figure 1: The DVCS (a) and the Bethe-Heitler (b and c) processes. 



The QCD interpretation of DVCS is based on the two diagrams shown in Fig. |^. In the 
presence of a hard scale, the DVCS scattering amplitude factorises [@,|3|,|^] into a hard scattering 
part calculable in perturbative QCD and parton distributions which contain the non-perturbative 
effects due to the proton structure. The DVCS process is similar to diffractive vector meson 
electro-production, a real photon replacing the final state vector meson. Recent measurements 
of diffractive vector meson production at HERA [[TT1-|T5|] indicate that Q^ is relevant for the hard 
scale of the interaction in QCD calculation. In comparison to vector meson production, DVCS 
avoids the theoretical complications and uncertainties associated with the unknown vector me- 
son wave function. However, even at Q"^ values above a few GeV^, non perturbative effects 









b) 



Figure 2: The two leading DVCS diagrams in a QCD picture. 



influence the predictions and have to be modeled. The DVCS cross section is suppressed rela- 
tive to that of vector meson production by the additional electromagnetic coupling. 

Calculations of the DVCS cross section have been published by Frankfurt, Freund and 
Strikman (FFS) \^ and by Donnachie and Dosch (DD) [[17p. Both contain "soft" (non- 
perturbative) and "hard" contributions. The soft part in the FFS prediction is based on the 
aligned jet model [p^, whereas reggeon and soft pomeron exchanges are considered in the 
DD prediction. The hard contribution in FFS is calculated in perturbative QCD. The hard 
contribution in DD is based on a dipole model where all parameters are determined from pp 
and 7*p cross section measurements. These predictions only provide the scattering amplitude 
at t = train — —m^Q^/W^, where |t| is the squared momentum transfer to the proton, nip 
the proton mass and W the invariant mass of the 7*^9 system. In both cases an exponential 
t-dependence, e^'''*', is assumed. 

For the DVCS process as well as for the vector meson production, the transition from a vir- 
tual photon to an on shell particle forces the fractional momenta of the two partons involved to 
be unequal and imposes a correlation. Therefore, the cross section calculation necessitates the 
use of skewed parton distributions [0-|^]. The difference in fractional momenta, the skewedness, 
becomes important at high Q^ values or high vector meson masses. In particular SPDs have 
been introduced in order to reconcile QCD calculations with the T diffractive photo-production 

i^. For vector meson production SPDs appear quadratically in 
A unique feature of DVCS is that they appear linearly in the 



measurements at HERA [ ]T9| , 
the cross section expression, 
interference term with the Bethe-Heitler process. Therefore the DVCS measurement, in con- 
trast to vector meson production, offers a particularly suitable channel to extract skewed parton 
distributions [|T^,p]- 



In this paper, the elastic DVCS cross section measurement at HERA is presented differen- 
tially in g2 and W in the Q^ range from 2 to 20 GeV^ W from 30 to 120 GeV and \t\ below 
1 GeV^. 



2 HI Detector, Kinematics and Monte Carlo Simulation 

A detailed description of the HI detector can be found in [E^. Here only the detector compo- 



nents relevant for the present analysis are shortly described. The SPACAL [Q, a lead - scin- 
tillating fibre calorimeter covers the backward^] region of the HI detector (153° < 6 < 177.5°). 
Its energy resolution for electromagnetic showers is a{E)/E ~ l.lVo/ ^jE/GeV © 1%. The 
uncertainty on the alignment of the calorimeter corresponds to an uncertainty of l.Smrad on 
the scattered positron polar angle. The liquid argon (LAr) calorimeter (4° < 6* < 154°) 
is situated inside a solenoidal magnet. The energy resolution for electromagnetic showers is 
a{E)/E ~ 11%/ y^E/GeV as obtained from test beam measurements p^. The major com- 
ponents of the central tracking detector are two 2 m long coaxial cylindrical drift chambers, the 
CJC, with wires parallel to the beam direction. The measurement of charged particle transverse 
momenta is performed in a magnetic field of 1.15 T, uniform over the full tracker volume. The 
forward components of the detector, used here to tag hadronic activity at high pseudo-rapidity 
(5^r]^7), are the forward muon spectrometer (FMD) and the proton remnant tagger (PRT). 
The FMD, designed to identify and measure the momentum of muons emitted in the forward 
direction, contains six active layers, each made of a pair of planes of drift cells. The three layers 
between the main calorimeter and the toroidal magnet can be reached by secondary particles 
arising from the interaction of particles scattered under small angles hitting the beam collima- 
tors or the beam pipe walls. Secondary particles or the scattered proton at large \t\ can also be 
detected by the PRT, located at 24 m from the interaction point and consisting of double layers 
of scintillator surrounding the beam pipe. The trigger used is based on the detection of a cluster 
in the electromagnetic section of the SpaCal calorimeter with an energy greater than 6 GeV. 

The data were obtained with the HI detector in the 1997 running period when the HERA 
collider was operated with 820 GeV protons and 27.5 GeV positrons. The sample corresponds 
to an integrated luminosity of 8 pb^^. 

The reconstruction method for the kinematic variables Q^ and s-Bjorken relies on the polar 
angle measurements of the final state electron, 9^, and photon, 6^^, (double angle method): 

q2 ^ ^^2 sin g^(l + cosge) 2) 

° sin 6j + sin 9^ — sin (6^ + 9^) 
Eq sin 9.^ + sin 9e + sin (6'e + 9^) 
Ep sin 9^ + sin 9e — sin {9e + 9^) 

W"^ = — (1-x) (4) 

X 

where Eq and Ep are the electron and the proton beam energies, respectively. In case a vertex 
cannot be reconstructed, the nominal position of the ep interactions is taken. The variable t, the 
square of the four-momentum transfer to the proton, is very well approximated by the negative 
square of the transverse momentum of the outgoing proton. The latter is computed as the vector 
sum of the transverse momenta of the final state photon pt^ and of the scattered positron pt^: 

t^-\Pu+Ptf. (5) 



' HI uses a right-handed coordinate system with the z axis taken along the beam direction, the +z or "forward" 
direction being that of the outgoing proton beam, the "transverse" directions are perpendicular to the z axis. The 
polar angle is defined with respect to the z axis. The pseudo-rapidity is defined by 77 = — In tan 9/2. 



Monte Carlo simulations are used to estimate the corrections to be applied to the data for the 
acceptance and resolutions of the detector. The generated events are passed through a detailed 
simulation of the HI detector and are subject to the same reconstruction and analysis chain as 
the data. The DVCS process is simulated (for more details see [^) according to the predicted 
cross section of FFS [ p^ , which includes the DVCS process, the Bethe-Heitler process and their 
interference^. Photon radiation from the incoming positron has been included in the simulation 
in the coUinear approximation. The Monte Carlo generator COMPTON 2.0 [ PD[ ] is used to 
simulate Bethe-Heitler events. To simulate the background sources (see section 3), diffractive 
p, uj and meson events are generated with the DIFFVM Monte Carlo [|3T|], dilepton production 
through a photon-photon interaction is simulated using the GRAPE program 



3 Event Selection 



The cross section of the Bethe-Heitler process, proceeding via Bremsstrahlung from the positron 
lines, is the largest when the positron and the photon are both produced in the backward direc- 
tion. In the DVCS case, the final state photon does not originate from the positron and therefore 
the ratio of DVCS over BH process is expected to increase when the photon is found in the 
forward direction. The analysis is thus restricted to the case where the photon is detected in 
the central or in the forward parts of the detector, i.e. in the LAr calorimeter. A data sam- 
ple, dominated by Bethe-Heitler events, is used as a reference sample to monitor the detector 
performance and its simulation. Two event samples are selected. 

• Enriched DVCS sample: The photon candidate is detected in the LAr calorimeter and 
the positron candidate in the SpaCal calorimeter. Both DVCS and Bethe-Heitler processes 
are expected to contribute to this sample with similar magnitudes. 

• Control sample: The photon candidate is detected in the SpaCal calorimeter and the 
positron candidate in the LAr calorimeter. The contribution of DVCS to this sample is 
negligible. 

The event selection is based on the detection of exactly two electromagnetic clusters, cor- 
responding to the final state photon and positron. One cluster is required to be detected in the 
SpaCal calorimeter with energy larger than 15 GeV and the other one in the LAr calorimeter 
(25° — 145°) with a transverse momentum pt > 1 GeV. Events with more than one track are 
rejected. Events with one track are only kept if the track is associated to one of the clusters 
which hence identifies the positron candidate. If no track is reconstructed, the SpaCal cluster 
is assumed to be the positron. In order to reject inelastic and proton dissociation events, no 
further cluster in the LAr calorimeter with energy above 0.5 GeV is allowed and the absence of 
activity above the noise level in forward detectors PRT and FMD is required. The influence of 



^ The free parameters have been set in the simulation to the following values: the t-slope parameter h — 
7GeV~^ the phase of the QCD ampHtude, -qQCD = 1 - f (0.176 + 0.33 log Q^) |||] and the sensitivity to 
the skewedness of the parton densities R-^ — 0.55 Jig], The proton structure function as extracted from the HI 
data 1129] has been used. 



QED radiative corrections is reduced by requirements on the longitudinal momentum balance^]. 
To enhance the DVCS signal with respect to the Bethe-Heitler contribution, and to maintain a 
large detector acceptance, the kinematic domain is explicitly restricted to: 2 < Q^ < 20 GeV^, 
|t| < 1 GeV^ and 30 < W < 120 GeV. It has to be noted that for the BH process, the Q^ and 
W variables cannot be associated with the photon virtuality and the 7*p center of mass energy, 
respectively. 



3.1 Control sample 

This sample of 338 events is dominated by the Bethe-Heitler process. Due to the large scatter- 
ing angle of the positron, the DVCS process is suppressed to negligible levels. In order to have 
a control of the detector response in the same energy and angular ranges as for the enriched 
DVCS sample, the kinematic cuts in Q"^ and W are applied to this sample, treating the photon 
candidate in SpaCal as the scattered positron and the positron candidate in the LAr calorime- 
ter as the photon. Background contributions from inelastic Bethe-Heitler events, diffractive p 
meson production and electron pair production have to be considered. The contribution of in- 
elastic Bethe-Heitler events is estimated to be 7.7 ± 3.8 %. The diffractive electro production of 
p mesons constitutes a background when one of the decay pions is mis-identified as a positron 
in the calorimeter and the other pion is undetected, while the positron scattered into the SpaCal 
calorimeter is taken to be the photon. The elastic production of electron pairs by photon-photon 
processes e^p — > e^e'e^p contributes to the background when only two of the three leptons are 
detected. These processes have been simulated using the Monte Carlo programs COMPTON 
2.0 (both for elastic and inelastic BH), DIFFVM and GRAPE respectively. In the left column of 
Fig. ^, event distributions of the control sample are shown. The data are compared to the sum 
of the absolute MC predictions which includes the Bethe-Heitler process, elastic p production 
and elastic dilepton production, all normalised to luminosity. A good description of the data 
by the sum of the different MC samples is achieved, showing that the detector response is well 
described by the simulation. 



3.2 Enriched DVCS sample 

This sample (172 events) is found to be dominated by the DVCS contribution although the 
contribution of the Bethe-Heitler process is non negligible. An important contamination to 
DVCS elastic candidates is due to the DVCS process with proton dissociation: 

e+ + p ^ e+ + 7 + F, (6) 

when the decay products of the baryonic system Y are not detected in the forward detectors. 
The sum of non elastic DVCS and BH contributions has been estimated to be 16 ± 8 % of the fi- 
nal sample p7|], based on the fraction of events with proton dissociation tagged by the forward 
detectors and the detection efficiency of the forward detectors for proton dissociation events 
using the DIFFVM Monte Carlo. The other sources of background to be considered are due 



^ The quantity ^ i? — P^ is required to be above 45 GeV. E denotes the energy and P^ is the momentum along 
the beam axis of the final state particles. The sum is calculated for the final state positron and photon. 
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Figure 3: Event distributions of the control sample (left) and of the enriched DVCS sample 
(right), a-b) energy of the cluster in the LAr calorimeter, c-d) polar angle of the cluster in 
the LAr calorimeter, e-f) coplanarity, i.e. difference of the azimuthal angle of the positron and 
photon candidates. The error bars on data points are statistical. Control sample: the cluster in 
the LAr calorimeter corresponds to the positron candidate. The data are compared to the sum 
of the predictions for the Bethe-Heitler process, elastic dilepton production and diffractive p 
production. All predictions are normalised to luminosity. Enriched DVCS sample: the cluster 
in the LAr calorimeter corresponds to the photon candidate. The data are compared to the sum 
of the predictions for the e^p -^ e'^'jp reaction according to EES, added to uj and diffractive 
backgrounds. The backgrounds and the BH contribution (shown on top of the backgrounds) are 
normalised to luminosity whereas the DVCS prediction is normalised in such a way that the 
sum of all contributions is equal to the total number of events. 



to diffractive uj and production, with decay modes to final states including photons (directly 
or from vr" decay) or K\. The main contributions originate from uj -^ 7r°7 and — > K^K^ 
followed by K^ -^ tt^tt". The background arising from tt" production, with the decay photons 
reconstructed in a single cluster, in low multiplicity DIS, is estimated from data. According to 
the MC simulation the two decay photons are recognised, in this energy range, in 20 % of the 
events in separate clusters. The selection has thus been extended to events with at least two 
clusters in the LAr calorimeter of energies exceeding 0.5 GeV. The two photon invariant mass 
spectrum is in agreement with expectation from diffractive oj and (p production. Background 
from Ti^ production is therefore inferred to be negligible. Figure ^ shows data distributions in 
comparison to the sum of the predictions according to the FFS calculation and the diffractive uj 
and backgrounds. The BH contribution in the FFS prediction and the uj and backgrounds 
are normalised to luminosity. The DVCS contribution in the FFS prediction is here normalised 
such that the sum of all contributions is equal to the total number of events in the data. The pure 
Bethe-Heitler contribution is also shown. The DVCS signal exhibits different kinematic distri- 
butions from the Bethe-Heitler contributions, in particular in the polar angle of the LAr cluster 
(Fig. ^-d) and in coplanarity (Fig. ^-f), which is defined as the difference of the azimuthal 
angles of the two clusters and is related to the p^-balance of the positron-photon system. The 
coplanarity distribution is found to be broader (rms = 5.3°) for the sum of the contribution in 
the enriched DVCS sample than for the Bethe-Heitler dominated control sample (rms = 3.2°). 
This is attributed to the electromagnetic nature of the BH process which implies a steeper t 
distribution than the DVCS signal. 



4 Cross section measurement and model comparison 

To extract the cross section, the data of the enriched DVCS sample have been corrected for 
detector effects, acceptance and for initial state radiation of real photons from the positron line 
using the Monte Carlo simulation. The bin size has been chosen according to the statistical 
accuracy and is large with respect to resolutions in Q^ (12 %) and W (6 %). The contamination 
of inelastic BH and DVCS events with proton dissociation (16 ± 8 %) is subtracted statistically. 
The background contributions from diffractive uj and production, estimated to be 3.5 % on 
average and below 6 % in all bins, have been subtracted bin by bin. 

The e'^p -^ e^'jp cross section is presented in Fig. ^ and in Table [1] differentially in Q^ and 
W. The data are compared to the estimate of the pure Bethe-Heitler contribution. The total cross 
section is dominated by the DVCS process at small W values and the Bethe-Heitler process at 
large W values. With the present precision, the Q^ slopes of the two processes appear similar. 
The limited resolution and statistics do not allow the cross section measurement differentially 
in t and the extraction of the slope from the present data. The main contribution, of 8 %, to the 
systematic error arises from detector effects and is due to the uncertainty on the measurement 
of the angle of the scattered positron. Other detector related errors are estimated to be around 
or below 2 %. The second largest systematic error arises from the estimate of the contamination 
of non elastic BH and DVCS (8 %). The total systematic error is found to be around 15 %. 

In the leading twist approximation the contribution of the interference term to the cross 
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Figure 4: Differential cross section measurements for the reaction e^p -^ e+7p as a function 
of Q^ (a) and W (b). The inner error bars are statistical and the full error bars include the 
systematic errors added in quadrature. The hatched histogram shows the contribution of the 
Bethe-Heitler process to the reaction, where, however, Q^ and W do not correspond to the 
photon virtuality and the YP center of mass energy, respectively. 

section is proportional to the cosine of the photon azimuthal angleQ. Since the present measure- 
ment is integrated over this angle, the overall contribution of the interference term is negligible. 
Therefore the Bethe-Heitler cross section can be subtracted from the total cross section in order 
to obtain the DVCS cross section. The DVCS e^p -^ e'^'jp cross section is then converted to a 
DVCS 'j*p -^ 7P cross section using the equivalent photon approximation (as in [[TT|]). 



g2 [GeV^; 


da^'^P- 


^^^pydQ^ [pb/GeV^" 


W [GeV 


da'^P- 


^^+pi/dW [pb/GeV 


[2.0,4.0] 


37.6 


±5.3 ±5.1 


[30, 60] 


1.06 


±0.15 ±0.10 


[4.0,6.5] 


8.0 


±1.3 ±1.1 


[60, 80] 


1.05 


±0.19 ±0.11 


[6.5,11.0] 


2.87 


±0.46 ±0.35 


[80, 100] 


1.58 


±0.23 ±0.16 


11.0,20.0 


0.61 


±0.12 ±0.11 


[100, 120] 


1.66 


±0.33 ±0.21 



Table 1: Differential cross sections for the reaction: e^p -^ e^p'j, as a function ofQ^ and W, 
in the kinematic domain 2 < Q^ < 20 GeV^ 30 < W < 120 GeV and \t\ < 1 GeVl The first 
errors are statistical, the second systematic. 

The 7*p cross section for the DVCS process is shown in Fig. |] and given in Table |^ as a 
function of Q"^ for H^ = 75 GeV, and as a function of W for Q"^ = 4.5 GeV^. The systematic 
errors on the 7*p cross section are due to the propagation of the systematic errors on the e~^p 
cross section combined with the bin center corrections error (7 %). The data are compared with 
the predictions by FFS and DD. The shape of the data is well described by both calculations 
both in Q^ and W. The absence of predictions for the t-slope leaves an uncertainty on the 



"^The photon azimuthal angle is defined as the angle between the plane formed by the incoming and scattered 
positrons and the 7* proton plane. 
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Figure 5: Cross section measurements for the ^*p — > 7^9 DVCS process as a function of Q^ 
(a) and W (b). The data are compared to the theoretical predictions of FFS [|T^ and DD [[T7|]. 
The band associated to each prediction corresponds to a variation of the assumed t-slope from 
5 GeV^^ (upper bound) to 9 GeV^^ (lower bound). The inner error bars are statistical and the 
full error bars include the systematic errors added in quadrature. 



normalisation of the theoretical models. The band associated to each prediction corresponds to 
a variation of the t-slope of 5 < 6 < 9 GeV^^, covering the measured range for light vector 
meson production [ [TT| , p^ ]. Both predictions are consistent with the data within this uncertainty. 
It is noted that these data provide constraints also for recent NLO calculations, which invoke 
skewed parton distributions []33|]. 



Q2 [GeV^; 


11.0 ±2.4 ±2.5 


W 


GeV 


(ji-p^jp [nb; 


3.0 




45 


4.33 ±0.64 ±0.54 


5.25 


3.8 ±1.1 ±1.0 




70 


5.51 ±1.34 ±0.86 


8.75 


2.43 ±0.66 ±0.54 




90 


8.9 ±2.2 ±1.7 


15.50 


0.64 ±0.30 ±0.28 




110 


4.8 ±4.0 ±2.6 



Table 2: Measured cross section for the elastic DVCS process YP ^ 'jp as a function ofQ^ 



for W = 75 GeV and as a function ofW for Q"^ 
errors are statistical, the second systematic. 



4.5 GeV^ , both for \t\ < 1 GeV^ The first 



5 Conclusion 

The DVCS process has been studied in the kinematic region 2 < Q^ < 20 GeV^, 30 < ly < 
120 GeV and |t| < 1 GeV^ using a data sample taken by the HI detector and corresponding to 
an integrated luminosity of 8pb^^ The cross section for the reaction e^p -^ e+7p has been 
measured for the first time and presented differentially in Q'^ and W. The DVCS process is 
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observed to dominate over the Bethe-Heitler process for W ^ 100 GeV. The 7*p DVCS cross 
section has been extracted and compared to the QCD based predictions [[T6|,p^] which both 
describe the measured Q'^ and W distributions within errors. 
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